Lie Symmetries and Exact Solutions 
of the Generalized Thin Film Equation 

Roman CHERNIHA % Phil BROADBRIDGE and Liliia MYRONIUK " 

Institute of Mathematics, NAS of Ukraine, 
3 Tereshchenkivs'ka Street, 01601 Kyiv, Ukraine 
E-mail: cherniha@imath.kiev.ua 
^ School of Engineering and Mathematical Sciences, 
La Trobe University, Bundoora VIC 3083, Australia 
E-mail: P. Broadbridge@latrobe .edu. au 
^ Faculty of Mathematics, 
Lesya Ukramka Volyn National University, 
13 Voly Avenue, 43025 Lutsk, Ukraine 
E-mail: Liliia_Myroniuk@univer.lutsk.ua 

Abstract 

A symmetry group classification for fourth-order reaction-diffusion equations, allowing for 
both second-order and fourth-order diffusion terms, is carried out. The fourth order equations 
are treated, firstly, as systems of second-order equations that bear some resemblance to systems 
of coupled reaction-diffusion equations with cross diffusion, secondly, as systems of a second- 
order equation and two first-order equations. The paper generalizes the results of Lie symmetry 
analysis derived earlier for particular cases of these equations. Various exact solutions are con- 
structed using Lie symmetry reductions of the reaction-diffusion systems to ordinary differential 
equations. The solutions include some unusual structures as well as the familiar types that reg- 
ularly occur in symmetry reductions, namely self-similar solutions, decelerating and decaying 
traveling waves, and steady states. 

1. Introduction 

We consider the fourth order nonlinear partial differential equation (PDE) of the form 

ut = + [D{u)u^]x + F{u), (1) 

where K, D and F are arbitrary smooth functions ( hereafter the subscripts t and x denote 
differentiation with respect to these variables). Eq. ([T]) generalizes a wide range of the known 
scalar reaction-diffusion equations arising in applications. The case with K identically zero 
and D{u) > for almost all u, is the case of second-order react ion- diffusion which has already 
been widely studied in many practical contexts including combustion, population dynamics, 
population genetics, neurobiology, biological cellular growth and adsorptive porous media, e.g. 
PJ-^. Hereinafter we assume that K{u) is not identically zero, so that the governing equation 
is of the fourth order, including a fourth-order diffusion term when K is non-negative. The 
simplest equation of the form ([1]), with F = and D = 0, follows from the approximations of 
lubrication theory to describe thin films of a Newtonian liquid dominated by surface tension 



effects. The thin film equations are an active area of research (see [S]^[H] and papers cited 
therein). The equation 



ut = -[u 




X 



(2) 



with the non-negative parameter 7 was introduced in [9]. 7 = 8 describes a classical thin film 
of Newtonian fluid, as reviewed in [10]. 7 = 1 occurs in the dynamics of a Hele-Shaw cell jll] 
and 7 = 2 arises in a study of wetting films with a free contact line between film and substrate 



One important generalization of Eq. ([2]), which is also a particular case of Eq. ([T]), reads as 



where /i is a positive parameter (arbitrary positive coefficients of each term can be set to one by 
rescaling variables). Eq. with 7 = can be considered as a semilinear limit of the classical 
Cahn-Hilliard model of phase separation [12], which is also widely studied (see [13] and the 
papers cited therein). The linear case 7 = /z = also follows from a small-slope approximation 
to metal surface evolution, with surface-diffusion and evaporation-condensation represented by 
fourth-order and second-order diffusion terms [Tll[T5|. The case 7 = /i = 3 arises in a study of 
capillary instability of axisymmetric thin films |16] . 

Several papers are devoted to the construction of exact solutions of the thin film equations 
02]) by Lie symmetry reductions and their generalizations [I7|-[2l], or by searching for special 
invariant finite vector spaces of solutions [25]. The symmetry classification is extended here to 
include a reaction term. In some circumstances, the reaction term F should naturally arise in 
fourth-order transport equations with a role similar to that in second-order react ion- diffusion. 
For example, a particular case of Eq. ([T]) with F{u) = u occurs as the limiting case of the 
unstable Cahn-Hilliard equation [26] [27] . In fabricated metal surface evolution, a positive 
source term may represent ion beam sputtering [28] and a negative source term may represent 
chemical decay or evaporation [29]. Other examples of equations of the form (JT]) arising in 
applications and having a reaction term, are presented in [27]. 

The first aim of this paper is to describe all possible Lie symmetries, which Eq. ([T|) can ad- 
mit depending on the function triplets [K, D, F), i.e. to solve the so-called group classification 
problem, which was formulated and solved for a class of non-linear heat equations in the pio- 
neering work in [30]. This problem for the second-order react ion- diffusion equation was solved 
in [21] (see also [52| [55] where the problem is solved for the general reaction-diffusion-convection 
equation). Note that the most general results concerning non-classical (Q-conditional) symme- 
tries of reaction-diffusion equations were obtained in [34] [36]. 

It should be noted that we shall not directly search for Lie symmetries of Eq. ([1]) but replace 
one scalar equation by an equivalent cross- diffusion system of equations. Using the symmetries 
found, we construct exact solutions of Eq. ([T]) with such triplets {K,D,F), which arise in 
applications and compare the results obtained with those derived earlier. 

Ovsiannikov's method of Lie symmetry classification of differential equations [37] is based 
on the classical Lie scheme and a set of equivalence transformations of a given equation. The 
formal application of this method to equations containing several arbitrary functions (Eq.([T]) 



Ut = - [u^U^xx]x + [U^U^] 



x-i 



(3) 
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contains three arbitrary functions) usually leads to a large number of equations admitting non- 
trivial Lie algebras of invariance. Our approach of Lie symmetry classification of differential 
equations is based on the classical Lie scheme and on finding and then making systematic use 
of the sets of local transformations that reduce any differential equation with a Lie algebra of 
invariance, to one given in the relevant list, that is representative of each equivalence class. 
This approach has earlier been applied also for reaction-diffusion systems |40]. 

The paper is organized as follows. Section 2 is devoted to a complete description of Lie 
symmetries of Eq. ([1]), i.e. all possible Lie symmetries, which this equation can admit depending 
on the form the functions K, D and F, are found. In most applications, K and D are non- 
negative so that the diffusive transport processes are dissipative [H]. However, in solving 
a group classification problem, it is more common to allow the coefficient functions to be 
arbitrary smooth functions. In Section 3, the symmetry reductions and some exact solutions 
are constructed for particular cases of Eq. ([1]) that are likely to be useful in applications. The 
main results of the paper are summarized in the last section. 

2. Lie symmetry of Eq. ([1]) 

Firstly, we note that Eq. ([T]) can be reduced to the system 

ut = -[K{u)v^J^.^ + [D{u)u^].^ + F{u), 

(4) 

= Uxx — V 

by the substitution v = Uxx, v = v{t,x). Since Ux is not relabeled as a separate variable, point 
symmetries of this system do not include contact symmetries of the original single equation 
([T]). Since second-order contact symmetries do not exist point symmetries of this system 
could include nothing more than extensions (prolongations) of point symmetries of the single 
fourth order equation for u{x,t). However, it is convenient to analyse system (jl]) as a cross- 
diffusion system, in which the second equation contains the time variable t as a parameter. The 
motivation follows from the well-known fact that application of Lie's algorithm in the case of 
high-order equations leads to very cumbersome formulae. Moreover, each Lie symmetry of Eq. 
([T]) can be easily established from one of system (jlj) and there is one-to-one correspondence 
between solutions of (jl]) and ([1]). Thus, we shall investigate system (jl]) instead of Eq. (jl]). 

Lie's classical algorithm has been implemented within several computer-algebra routines. 
The classification problem may be solved by iterating a symmetry-finding program such as 
the REDUCE-based program DIMSYM |43j or MAPLE-based program DESOLV [11] that flag 
special cases of the coefficient functions that may give rise to additional symmetries. The 
REDUCE-based program CRACK ^45j can fully solve some classification problems automati- 
cally. However, these programs may produce many special cases of equations with additional 
symmetries, which are equivalent up to the correctly specified local substitutions. Finding and 
systematically using the sets of such substitutions often leads to a significant practical reduc- 
tion in the number of special cases that admit additional invariance. Thus, the problem under 
investigation here will be tractable without the assistance of a computer. 

Now we formulate the main theorem, which presents the classification of special forms of 
system (jl]) (with K not identically 0) having additional symmetries. 
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Theorem 1 All possible maximal algebras of invariance (up to equivalent representations gen- 
erated by transformations of the form ^) of system ^ for any fixed triplet {K, D, F) are 
presented in Table 1. Any other system of the form with non-trivial Lie symmetry is re- 
duced by a local substitution of the form 

t-^Cot + Cie^^\ 
X — i- Csx, 

(5) 

u^Ci + C^t + C^e^^'u, 



to one of those given in Table 1 (the constants C with subscripts are determined by the form of 
the system in question, some of them necessarily being zero in all particular cases). 

Remark 1. All the systems listed in Table 1 are inequivalent up to arbitrary local substitution 
(not only of the form ([5])!). This can be shown using the same approach as that used in [32] 
for the general reaction-diffusion-convection equation. 
Proof. 

According to the classical Lie scheme |1S]^[1H], we consider system (jl]) as the manifold 
(•S*!, 5*2) determined by the restrictions: 

Si = -ut - [K{u)v,], + [D{u)u,], + F{u) = 0, 

(6) 

5*2 = Uxx — f = 

in the space of the variables t, x, u, f , Ut, Vt, Ux, Vx, Uxx, Vxx- The maximal algebra of invariance 
(MAI) of this system is generated by infinitesimal operators of the form 

X = ^°(t, X, u, v)dt + i^{t, X, u, v)dx + ri^{t, x, u, v)du + rj^{t, x, u, v)d^, (7) 

where the functions C,^, rj^, rf are to be determined. In order to determine these unknown 
functions one needs to use the invariance conditions 



X^S, = X\-Ut - [K{u)vx]x + [D{u)ux]x + F{u)) 



= 0, 

51=0,52=0 



X^S2 = X\uxx-v) 



= 0, 

51=0,52=0 



(8) 



where is the second prolongation of the operator X: 

X^ = ^"(t, X, u, v)dt + S}{t, X, u, v)dx + ffit, X, u, v)du + Tf'it, x, u, v)dv-\- 

+ r^idu^ + rfodvt + r^A^ + (9) 
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Table 1. Lie symmetries of 03]) 



Case 


Systems of the form (jl]) 


Restrictions 


Basic operators of MAI 


1. 


Uxx — V = 




Pt, Px) 

Di = Atdt + xdx - 2vdy 


2. 


Uxx -V = Q 


+ A V 


Pt, Px, 

D2 = 2(7 - 2fi)tdt + (7 - f-i)xdx 
+2idu - (7 - f^)vdu) 


3. 


Ut = -[u^v^]^ + c/K^J, + Am^/^-^+i 

Uxx — V = Q 


7' + /iV0 
+ A V 


Pt, Px, 

D3 = 2(7 - 2fi)tdt + (7 - fi)xdx 
+2{udu + (/i - 7 + l)vd^) 


4. 


Ut = —Vxx + duxx + Am In M 

Uxx — W = 


A ^ 


Pt, Px, 

Qi = e^\udu + vd^) 


5. 


Ut = -[e'^"t^x]x 

Uxx — f = 


77^0 


Pt, Px, 

Di, D2 with /i = 7 


6. 


Ut = -[u^v^]^ 

Uxx — V = 


77^0 


Pt, Px, 
Di, D3 with /i = 7 


7. 


Ut Uxx ~l~ dUxx 
Uxx - V = 


d^O 


Pt, Px, I = udu + vdy, 

X^ = P{t,x)d^ + Pxx{t,x)d,, 

where Pt + Pxxxx - dPxx = 


8. 


Ut = —Uxx 
Uxx - U = 




Pt, Px, I, Di, 

X~ = P{t,x)du + Pxx{t,x)d,, 
where Pt + Pxxxx = 



The coefficients with relevant subscripts Tt, rl , rfg, Tq^^, r2o, T^j^i, Tq2) ^loi '^11 i '^02 ^'^^ 
calculated by the well-known prolongation formulae (see, e.g. |l6]-|48j). 

Substituting operator ([9]) into system (|8]) and carrying out the relevant calculations, we 
obtain so called system of determining equations for finding the functions 77^, rf. This 

is an overdetermined system of PDEs that can be written in the explicit form: 

^x ~ ~ ~ ^t ~ ~ ~ ^xx ~ (lO) 
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vl = vlu = vlu = o (11) 
= vL + ivl - '^Qy (12) 

Kiu)i^?-Ae.) + K'iuW = (13) 
Dium-2Q + D'iuW = (14) 

Vl + Kiu)r^l,, + F{u){r,l - e°) - D{u)r,l^ - F\uW = (15) 

Of course, Eqs. f|T0|) - f|T2|) are linear and don't depend on the functions K, D and F, hence 
their general solution can be easily constructed: 

r]^ = R{t)u + P{t,x), 

T]^ = P^^{t,x) + iRit)-2a)v, 

where P{t,x), R{t) are arbitrary smooth functions, a and Xq are arbitrary constants. Eqs. 
(fT3|) - (fT5l) form the system of classification equations. Its general solution under assumption of 
arbitrarily given functions K, D and F generates the invariance algebra that consists only of 
generators of translations in x and t with the basis 

The algebra with this basis is called the trivial Lie algebra of the system (jlj) (note that other 
authors, instead use 'kernel of the basic Lie groups' [37] or 'the principal Lie algebra' p4] in this 
context). Thus, we aim to find all triplets [K, D, F) that lead to extensions of the trivial Lie 
algebra generated by ( ITTll . This means that one needs to solve Eqs. (IT3l)-(IT5l) with coefficients 
f lT6|) . The crucial step in solving this task is to analyze differential consequences of Eqs. flT3|) - 
f lT^ with respect to the variable x. Since these consequences take the form K'{u)Px{t,x) = 
and D'{u)Px{t,x) = 0, respectively, one arrives at two basic cases 

i [K'{u)]^ + [D'{u)]^ ^0 ii K'{u) = D'{u) = 0. 

Consider case i. In this case Px{t, x) = 0, so that Eqs. f|T3l) - f|T5l) take the form 

K{u)[^tit) - 4a] + K'{u)[R{t)u + P{t)] = 0, (18) 

D{u) [e°(t) - 2a] + D'{u) [R{t)u + P{t)] = 0, (19) 
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Rt{t)u + Pt{t) + F{u)[R{t) - ^°(t)] - F'{u)[R{t)u + P{t)] = 0. (20) 

Setting R{t) = P{t) = 0, we immediately arrive at case 1 of Table 1. In fact, Eqs. (1T8|) - (!20I) 
with R{t) = P{t) = and non-zero K{u) are equivalent to 

e?(t) = 4« ^ 0, D{u) = 0, F{u) = 

and K{u) is an arbitrary smooth function. This means that the triplet {K{u), 0,0) forms the 
system from case 1 of Table 1 and the coordinates of the infinitesimal operator ([7]) take the 
form 

= Aat + tr), f^=aX + Xn, 
Tj^ =7]^ = 0, 

where a, to and xq are arbitrary parameters. Operator ([7]) with coordinates f l2T|) generates 
exactly the operators Pt (for = l,a = = 0), Px (for = a = 0,xq = 1) and Di (for 
a = 1, Xo = to = 0) listed in case 1 of Table 1. 

\iR^{t) + P'^{t) ^ then two possible subcases arise: ia. R{t) 7^ and ib. R{t) = 0, P{t) ^ 

0. 

Consider subcase ia. Integrating Eq. ( ITSjl as an ODE on the function K{u), one obtains 

m 

R{t) 



K{u) = k u+^-^ , (22) 



where is a nonzero constant, which can be reduced to = 1 (by scaling time t — )■ kt) without 
losing generality. Since the function K must depend only on u, the restrictions 

et{t)=Aa-^R{t), 
Pit) = joR{t) 

are obtained. Here 7 and 70 are arbitrary constants. Thus, solving Eq. f[T5]) . we arrive at the 
power function 

K{u) = {u + ^o)^ (24) 
and restrictions (!23|) . 

To solve Eqs. ( !T9|) -( |20|) . we need to analyze two subcases ial D{u) 7^ and ia2 D{u) = 0. 
In subcase ial, Eq. (ITQll can be solved as an ODE on the function D{u), so that one obtains 

D(m) = c/(m + 7o)^ (25) 

and the condition 

e?(t) = 2a-/ii?(t), (26) 
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where (i 7^ and /i are arbitrary constants. 

If 7 7^ /i then solving Eqs. and and substituting the found functions R(t),P(t) 

and ^'^ into Eq. fl20|) . we arrive at the ODE 



{u + i^)F\u) + (7 - 2/i - l)F{u) = 0. (27) 

The general solution of Eq. f l27|) and the functions i^'(u) and -D(u) found above form the system 

ut = - [{u + 7o)^^^Jx + d[{u + 7o)^M J, + A(u + 70)"^"^+' , , 

(28) 

Uxx-V = 0, 

which is reduced to the system listed in case 3 of Table 1 by renaming m + 70 — )■ u. Substitut- 
ing the functions R{t),P{t) and into f lT6|) . we obtain the infinitesimal operator ([7]), which 
generates three basic operators listed in case 3 of Table 1. 

If 7 = /X then Eqs. fl25]) and fl2^ are compatible only under restriction a = 0. It turns 
out that this restriction doesn't lead to any new cases but to case 3 of Table 1 with 7 = /x. In 
fact, Eqs. ( !T8|) and (fT9|) have identical structure, hence D{u) = d{u + 70)''', d = const. The 
general solution of (!20|) has the form F{u) = X{u + 70)^^'*^ + Ai(m + 70), where A and Ai 7^ are 
arbitrary constants. The relevant coordinates of infinitesimal operator ([7]) take the form 

(29) 

^1 ^ ^_|^g -^i'>'*(m + 7o), = rie , 
where to? ''"1 ^^^d Xq are arbitrary parameters. Thus, the system 

Mt = -[(n + 7o)'^t;Jx + c?[(m + 7o)^Mx]x' + A(m + 70)^+'^ + Ai(m + 7o) 

(oUj 

Uxx — = 

is invariant under three-dimensional MAI with the basic operators 
P P 

g* = e-^^i*(9t + Ai((m + 7o)5. + t;9,)). ^^'^ 
However, we found the local substitution 

A17 ' 

" f ' N w (32) 
ti* = ve~^^\ 

which reduces system ( l30i) and Lie algebra ( l3Ti) to the system and Lie algebra listed in case 3 
of Table 1 with 7 = /i. 

The analysis of subcase ia2 is straightforward because Eq. (fT9|) vanishes for D{u) = while 
Eq. (|20l) can be treated in a similar way. In conclusion, we found that the system 

= -[(u + 7o)'^i;J^. + Ai(u + 7o), 

(33) 

Mxx - = 0, 
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is invariant under four-dimensional MAI with the basic operators 

= X\ ^{u + 7o)5„ + - 2)vd,. (^^^ 

Direct checking shows that system fl55]) and Lie algebra flM|) are reduced to the system and Lie 
algebra listed in case 6 of Table 1 if one applies substitution (!32|) . 

Examination of subcase ib is much simpler because Eqs. (fT8|) - (!20|) with R{t) = can be 
easily integrated. Finally, one obtains cases 2 and 5 of Table 1. 

To complete the proof we need to examine case ii K'{u) = D'iu) = 0, i.e. K{u) = k = const, 
D{u) = d = const. Since K{u) ^ 0, we can again set = 1 without losing generality. 

Thus, the classification equations f|T3|) - f|T4|) with coefficients f|T6|) can be essentially simplified 
and one obtains 

^O(t) = 4at + to, arf = 0. (35) 
The third classification equation f|T5|) takes the form 
Rt{t)u^ Pt{t,x) \ P^^^^{t,x) - dP^^{t,x)\ 

(36) 

+F(m) [i?(t) - 4a] - F'(m) [i?(t)M + P(t, x)] = 0. 

Differentiating Eq. fl36p with respect to x and m, we find the condition F"{u)Px{t,x) = 0. 
Hence two different subcases should be examined: 

iia F" = iib F" 0, P^{t, x) = 0. 

Consider subcase iia. Since F"{u) = we immediately obtain F{u) = Xiu + Aq so that the 
triplet of functions {K, D, F) is known. 

Substituting the function F into Eq. fl5B]) and splitting the obtained expression into two 
equations (with the variable u and without it) we arrive at R{t) = 4a\it + ri and the linear 
PDE 

Pt{t, x) + P^^^^it, x) - dP^^it, x) - AiP(t, x) + Ao(4aAit + n - 4a) = (37) 
to find the function P{t,x). Thus, the coordinates of infinitesimal operator (jTj) take the form 
f ° = 4at + to, = ax + xn, 

(38) 

ryi = (4aAit + ri)u + P{t, x),7f = {AaXit + ri - 2a)v + Pxx{t, x), 

where P{t, x) is the general solution of Eq. fl57|) . 

The last step is to take into account the second condition from 0351) . If c? 7^ then a = 
and, applying the relevant simplifications, we arrive at the case 7 of Table 1. 

If (i = then Eq. (137|) takes the form 

Pt{t, x) + P,,,,(t, x) - \iP{t, x) + Ao(4aAit + n - 4a). (39) 
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The corresponding system is 

ut = -Vxx + \iu + Ao, ^ ^ 

(40) 

Uxx — V = 

It turns out that system (l39l) is reduced to the form hsted in case 8 of Table 1 if one apphes 
the local substitutions 

u* = u - Xot, Ai = 

(41) 

V = V 



and 



u* = e-'^\u + j:^), Ai^O 

* Alt ^ ^ 

V = e If. 

Simultaneously, operator ([7]) with coordinates fl55]) is transformed in such a way that the basic 
operators Pt, Px, I, -Di and = P(t,x)du + Pxx(t,x)djj, hsted in case 8 of Table 1, can be 
easily derived. 

Consider subcase lib. Since Px{t,x) = 0, Eq. (136|) takes the form 

Rt{t)u + Pt{t) + F{u)[R{t) - 4a] - F'{u)[R{t)u + P{t)] = 0. (43) 
Differentiating Eq. fH5]) with respect to the variables t and u, we find the equation 

F"{u)[Rt{t)u + Pt{t)] = Ru{t). (44) 

Because Eq. has a simple structure we prefer to solve this equation and check when the 
solution obtained will satisfy Eq. (B3|) . We note that the special case Rt(t) = Pt(t) = doesn't 
lead to new results, so that Rf{t) + Pt{t) ^ 0. Moreover, since the function F depends only on 
u, the relations 

Pt{t)=^Rt{t), 

(45) 

Ru{t) = XRtit), 

where 7 and A 7^ are arbitrary constants, should take place. The general solution of Eq. 
with the coefficients fH5l) is 



F(m) = A(M + 7)ln(u + 7) + AiM + Ao, (46) 

where Aq and Ai are arbitrary constants. Now we substitute f H6|) and the general solution of 
the linear ODEs system ( H5|) into Eq. (143|) and find conditions when the obtained expression 
can be fulfilled. The simple calculations give 

R(t) = ne^*, P{t) = 7rie^* 
a = 0, Ao = A17. 
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So the system 



(48) 

Uxx — V = 

admits MAI generated by operator ([7]) with coordinates 

(49) 

f]^ = rie^^{u + 7), Iff = Tie^^v. 
Finally, the system (HHIl and operator ([7]) with (H9!) are simplified by the substitution 

u* = e^iu + 7), 

^ (50) 

f * = e A ti, 

so that the system and MAI listed in case 4 of Table 1 are obtained. 

Thus, the system of determining equations ( fTOl) — ( |T5l) is completely solved and eight dif- 
ferent systems of the form (jlj) have been found, which admit three- and higher- dimensional 
Lie algebras. Simultaneously, we have shown that all other systems admitting non-trivial Lie 
algebra are reduced to those listed in Table 1 by the substitutions of the form ( 132|) . pij) . ( 142|) 
and ( l50l) . One notes that all of these substitutions can be united to the form ([5]). 

The proof is now completed. ■ 

One easily notes that cases 2 and 3 of Table 1 generalize the results of Lie symmetry analysis 
for the Cahn-Hilliard equation derived in [191 El]- For example, the systems 

(51) 

Uxx — V = 

and 

Ut = -Uxx + d\u^u^^, 

(52) 

Uxx-U = 0, 

which are the particular cases of the corresponding systems from Table 1, admit the Lie sym- 
metry operators 

4/it5f -I- [ixd^ — 2((9„ -I- iivd^ (53) 

and 

4/it9t + /ia;9^-2(u(9„+ (/i + l)t;(9^) (54) 



respectively. The table also includes the results obtained in [21] for the equation Eq. ([T]) with 
K(u) = const and F{u) = 0. 
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To finish the Lie symmetry description we note that Eq. ([T]) can be reduced to an equivalent 
system of PDEs in different ways. One sees that the system 

ut = -[K{u)v^]^, + [D{u)u^]^ + F{u), 

= Ux — w, (55) 

= Wx — V 

by introducing new unknown functions v = v{t,x),w = w{t,x) can be obtained from Eq. ([T]). 
Since the system includes the first-order variable w = Ux, point symmetries of this system 
include contact symmetries of the original single fourth-order equation. System fl55l) is nothing 
else but a cross-diffusion system, in which the second and third equations contain the time 
variable t as a parameter. Thus, we may investigate also system fl55l) instead of Eq. ([1]). There 
is an essential difference between systems (jl]) and ([S5D because the second system contains the 
equations of different orders. In fact, according to the classical Lie scheme, MAI of this system 
is generated by the infinitesimal operator 

X = X, u, V, w)dt + i^{t, X, u, V, w)dx + y]^{t, x, u, v, w)du+ 
+ri^(t,x,u,v,w)dy + +'r]^(t,x,u,v,w)diu, 

where the functions rj^, rf,rf are to be determined. Applying the second prolongation 

of the operator fl56|) to system fl55l) and using the invariance conditions one can derive the 
determining equations for finding the functions r/^, rf^rf. It should be stressed that 

the relevant invariance conditions must take into account the differential consequences (with 
respect to the variable x) of the second and third equations of system (1^ . We omit cumbersome 



calculations and present the final result in the explicit form: 

— — it — iu — iv — ixx — (57) 

vl = vl = vL = vL = o (58) 

= vL + ivl - '^Qv (59) 

V' = Vl + ivl - (60) 

Kiu)i^^-AQ + K'iuW = (61) 

Diu)i^^-2ex) + D'iuW = (62) 

Vl + Kiu)vL.. + nu)K - it) - D{u)ri, - F'iuW = (63) 



Now we note that the determining equations obtained (of course, without Eq. (HI]) ) are 
equivalent to those for the system (jl]) so that no new Lie point symmetries or contact symmetries 
can be found. 
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3. Symmetry reduction and exact solutions 

Some of the systems presented in Table 1 are equivalent to known fourth order PDEs 
arising in applications. For example, the system listed in case 6 is nothing else but the thin 
film equation ([2]). Since the motivation to this study is to consider this type of equation with 
non-zero reaction terms, which naturally arise in some processes [26]- [29], henceforth we restrict 
our attention mainly to case 3 of Table 1. 

It is well-known that a Lie symmetry allows one to reduce the given PDE (system of PDEs) 
to an equation (system of equations) of lower dimensionality Here we reduce systems arising 
in case 3 of Table 1 to systems of ordinary differential equations (ODEs), furthermore these 
ODE systems are solved in particular cases and exact solutions of the initial PDE systems are 
constructed. Finally, these solutions are compared with those obtained by other authors. 

Consider the system arising in the case 3 of Table 1: 

(64) 

Uxx-v = 0, 

where 7^ + /i^ 7^ 0. The most general Lie symmetry operator of system flMl) has the form 

X = aiPt + a2Px + as-Ds = , , 

(65) 

= [ai + 2^3(7 - 2^)t]dt + [a2 + 03(7 - ^^)A^x + [2a2,u]du + [203(1 - 7 + ^)v\dv, 

where i = 1, 2, 3 are arbitrary constants. To construct the relevant ansatz one needs to solve 
the Pfaffian system of characteristic equations 

dt dx du dv 



— . (66) 
Oil -|- 2«3(7 — 2yu)i(: 0:2 + ^3(7 — yu)a; 2a-iU 203(1 — 7 + //)v 

The general solution of fl66l) essentially depends on the parameters ai, 0^2, 03, 7 and and five 
different cases occur. 

Case 1 03 = leads to the plane wave solutions of the form 

uj = aix — a2t, u = 4>{uj), v = ip{uj), (67) 

where and ip are new unknown functions. These functions should satisfy the ODE system 

-a2(f)' = -al(j)^-^[-f(j)'ilj' + (pip"] + dal(p>'-^[fi{(p')^ + cpcp"] + A^^-^+^z^, 

iP = a\<P". ^^^^ 

It should be noted that this system is equivalent to the fourth-order ODE 

+ ^a\<p-^<p'<p"' - da\<p^<p" - da\y.<p^~^ {uj) {<p'f - a2<p' - A^^-^+^m _ g. (gg) 

This equation is not integrable because there are no general solutions for this equation in terms 
of elementary functions and known special functions [49]. However, it can be noted that the 
special case with 7 = 3yU possesses the particular solution 

0(0;) = auj'^ , (70) 
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where a is a solution of algebraic equation 

afil - 2/i)a^'^ - dajf/a'^'' - aifx^af" - = 0. (71) 

Hence the system 



Uxx - V = 
possesses the exact solution 

u = a{aix — ^2^)'' , 



(72) 



(73) 



where a satisfies f l7T|) . 

Case 2 tta 7^ 0, 7 = 2yU 7^ 0, ai = leads to the ansatz 

— — — 2 

u = t, u = (j)(t)xt', v = il){t)xf^ . (74) 
The corresponding ODE system takes the form 

0' = -2(i - 1)(1 + f )02M^ + rfj( J + 1)0^^+^ + A0, 



(75) 



and can be rewritten as the single ODE: 

0' = --(- - 1)(^ - 1)0'^+' + d-{- + 1)0^+^ + A0. (76) 

It should be noted that this ODE with A = coincides with one derived in the recently 
published paper [24] for the fourth-order PDE, which is equivalent to f l64p with A = 0. However, 
system fl64|) with A 7^ cannot be reduced to one with A = so that the solutions presented 
below cannot be obtained from |24j . 

If /i = 1 or /i = 2 then the known solutions of the reaction-diffusion equations 

Ut = d[uux]x + Ati = 1, 

and 

Ut = d[u'^Ux]x + Am /i = 2 
are obtained because Uxxx = (see (ff4l)). If /i = — 2 then 

0(t) = Ce^' (77) 
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and there follows an exact solution 



u 

V 



X ' 



(7J 



of the system 



ut = -[:^v^]^ + d[:^u^]^ + Xu, 
Uxx -v = Q. 



(79) 



If /i 7^ 1; ±2 then two subcases, A = and A 7^ 0, should be separately considered. Both of 
them lead to the function (f){t) in the implicit form and one can be found from the transcendental 
equation 



In 



1 + 



d 



2(l-i)(^-l) 



d 



1 



it -to) 



(80) 



if A = and from the equation 



In 



(t) 



r{t) 



— a 







dz 



2A/i(t-to), 



(81) 



a 



if A ^ 0. 



2(1-m)(m-2)' 



(3 



4(1-m)(m2-4) 



Thus, the system 

Ut = -[u^^v^]^ + d[u''u^], 

Uxx — V = 

possesses the exact solution 



(82) 



(83) 



where (j){t) satisfies the equation (IHOj) . Note that the function (j){t) tends to if t — )■ 00 and this 
function blows up if t — )■ to- These properties follow from the simple analysis of (IHOj) . 
Analogously, the system 



Ut = -[u'^'^v^]^ + dlu^'u,,]^ + Xu, 
Uxx -V = 

possesses the exact solution (!83|l with 0(t) satisfying the equation (IHTl) . 



(84) 
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Case 3 as 7^ 0, 7 = 2/i 7^ 0, ai 7^ leads to the ansatz 

a; = °i , m = (j)[uj)e , f = '?/^(a;je , 

which reduces the initial system to the ODE system 

-a3fiuj(f)' + 2a30 = -ai02M-i[2/x0V' + + dai(f)''-^[fi{(f)')^ + (j)(f)"] + aiA0, 
^ = 0". 

System fl86|) is equivalent to the 4-th order equation 

+ 2/iai02/^-VV" - dai(f)''(l)" - dal^I(f)^''\uJ){(f)'f- 
—asfiojcf)' — (aiA — 20:3)0 = 0, 

which possesses the particular solution 



2 



Here, a must be a solution of the algebraic equation 
4(4 - _ ^)tt2/. _ 2rf;^2(^;^ ^ 2)a'' - A/i^ = 0, 

which is simply a quadratic equation for a'^. 
Thus, the cross- diffusion system 

ut = -[u^'^v^]^ + dlu^'uj^^ + Am, 
Uxx -V = 

has the stationary solution 



2 

u = axi^ 



5 



where a satisfies fl89l) . 

Case 4a3 7^0, 7 = /i7^0 leads to the ansatz 

02 _i _j_ 

a; = a; -I Int, u = 4>(u)t f", v = ip{u)t 1^ . 

The corresponding ODE system takes the form 

020' - 2^30 = -2a3/i0'^-i[/i0V' + H"] + 2a^iid(t)t'-^[ii{<p'f + 00"] + 2a3/iA0'^+\ 
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and is equivalent to the 4-th order equation 



(94) 

+^20' - 2a3/iA0^+i - 2«30 = 0. 

For ^2 7^ 0, the solutions u{x,t) are travelling waves whose speed decreases in proportion to 
and whose amplitude decreases in proportion to t~^/^. 

Equation ( 19^ is not integrable but we were able to find the particular solutions if = 1, 
a2 = 0: 

(P{u) = - A + _2_ sin(^^ + ^^)^ oi + de^-^ = (95) 

and 

</)M = -A + Cie^- + g^e-^ ^4-rf^2_A^o. (96) 



Thus, using ansatz ( ]92l) . we arrive at the solutions 

(97) 



^-L^ sin(0x+6»o) 



V 



t 

3^02sm(6»x+6lo) 



which is a spatial sinusoid for which amplitude varies in proportion to and 

' (98) 



9Ci a2 

U = : — 



of the system 

ut = -[uv^]^ + d[uu^]^ + Au^ 

(99) 

Uxx — V = 0. 

In formulae (1971) and (1951) . we can also shift the time t to t — or t + with the positive 
parameter to. The first shift leads to a solution having blow-up at t = to; the second leads to 
a solution that avoids singularity at t = and tends to as t approaches oo. Blow-up and 
extinction are interesting phenomena in some applications. 

Finally, case 5 03 7^ 0, 77^ 2/i, 7 7^ leads to the similarity reduction 

fJ. — 'y 1 1 — 7+^^ 

oj = xt2(7-2M) ^ u = (j){uj)t'y-^'' , V = ii{u)t -i-'^y- (100) 

and to the ODE system 

(/i - 7)w0' + 20 = -2(7 - 2/i)0^-i[70V' + 4>V\ + 2(7 - 2/x)#'^-i[/i(0')' + 00"] + 
+2(7 - 2/i)A0i-^+2^, (101) 
= 0". 
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The equivalent 4-th order equation has the form 



2(7 - 2/i)0T0^^ + 2(7 - 2^)70T-V>"' - 2(7 - 2^)rf0^0" - 2rf(7 - 2ij,)fi(f)''-\(f)'f+ 
+ - 7)a;0' - 2(7 - 2/i)A0i-^+2A. + 20 = 0. 

whose solutions are self-similar by a scaling invariance. Although this equation is again not 
integrable, one may try to construct particular solutions in the form of a high-order polynomial. 
For example, setting d = 0, 7 = 1, /i = 0, the exact solution 

(f)(uj) = —00^ + -A (103) 
^ 120 6 ^ ' 

is obtained. Thus, we arrive at the solution 



7/ — -J—^ -I- 5 \y- 
" 120 t 6^''' 

~ ToT 
of the system 

Ut = -[uvx]x + A, 

Uxx - V = 0. 

Remark 2. Using the ad hoc ansatz 

u = Mi) + M't)x + M't)x'^ + + 04(^)a;^ 

V = 202(t) + 6(p-i{t)x + 1204(^)a;^ 

solution fll04p can be generalized to the form 

u = ^ + 30^ -I- 900^ + 6750^ + |At+ 



(104) 



(105) 



(106) 



= 2( ^ + 45^ + + 



t y ' t ' 10 t ' 
where Cj {i = 0, 3) are arbitrary constants. 

Similarly, setting d = 0, 7 = 1, 1^ = \y the solution 

(f){io) = (^=00^ + —a)' (108) 
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can be derived. Thus, the system 

Ut = -[uVx]x + >^Vu, 

(109) 

Uxx — V = 
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possesses the exact solution 




(110) 



1 , 
V = —- h 

^ 10 i ^ 




We note that solutions f ll04p and (11 101) have been earlier obtained in [25j via the method 
of invariant subspaces. Indeed, the formulas (3.29) and (3.76) |25j contain (11041) and (IllOp . 
respectively. 



We have carried out a symmetry group classification for fourth-order reaction-diffusion 
equations, allowing for both second-order and fourth-order diffusion terms. The fourth order 
equation has been treated, firstly, as a system of second-order equations that bears some resem- 
blance to a system of coupled reaction-diffusion equations with cross diffusion, secondly, as a 
system of a second-order equation and two first-order equations. It turns out that both systems 
lead to the same result of symmetry group classification. Our paper generalizes the results of 
Lie symmetry analysis derived earlier for particular cases of Eq. ([1]) . Moreover, we were able 
to construct all possible Lie symmetries, which Eq. ([T]) can admit depending on the function 
triplets {K,D,F). This distinguishes our investigation from those that have focussed on Lie 
symmetry of particular cases of the given fourth-order evolution equation. On the other hand, 
our result is analogous to that derived in [31] and [32] for second-order evolution equations. 

To the best of our knowledge, there is only the recently published book [25], where exact 
solutions have been found for some equations of the form ([1]) with F ^ 0. Thus, a fourth-order 
nonlinear equation with the non-zero reaction term in the form of system (1M|1 was examined by 
applying the Lie symmetry reduction. Where possible, we have constructed exact solutions to 
the ordinary differential equations that were obtained from this reaction-diffusion system. The 
solutions include some unusual structures as well as the familiar types that regularly occur in 
symmetry reductions, namely self-similar solutions, decelerating and decaying traveling waves, 
and steady states. Many of the functional relationships between the two symmetry invariants 
are quite simple, involving polynomials, algebraic functions, logarithms, exponentials and sinu- 
soids. However, there are some that have been reduced only to the solutions of transcendental 
equations (see formulas (IHO!) and (IHT!) ). 

Finally, it should be also noted that the nonlinear fourth order ODEs obtained in section 
3 can be solved by numerical methods and therefore solutions of the relevant generalized thin 
film equations will be constructed. 
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